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In this letter we derive an action principle for Coulomb collisions in plasmas. Although no natural 
Lagrangian exists for the Landau-Fokker-Planck equation, an Eulerian variational formulation is 
found considering the system of partial differential equations that couple the distribution function 
and the Rosenbluth potentials. Exact conservation laws are derived after generalizing the energy- 
momentum stress tensor for second order Lagrangians and, in the case of a test-particle population 
in a given plasma background, the action principle is shown to correspond to the Langevin equation 
for individual particles. Being suitable for discretization, the presented action allows construction 
of variational integrators. Numerical implementation is left for a future study. 


Introduction - The principle of stationary action and 
the Noether theorem provide powerful tools for analyz¬ 
ing mechanics of particles and fields. Using them, it is 
possible to find a conservation law for any symmetry in 
a Lagrangian. Both tools also apply to discrete mechan¬ 
ics, allowing one to construct variational integrators (see 
e.g. [1, 2]) which typically conserve the discrete coun¬ 
terparts of the conserved continuous quantities. It is no 
wonder that physicists have had a long-lasting love affair 
with action principles. 

In plasma physics, simulations are entering an era 
where conservation of energy and momentum is not only 
a desired feature but a necessity: it would be difficult 
to trust the results of, for example, global turbulence 
simulations if numerical methods resulted in artificial 
damping or growth of turbulence. Superior numerical 
properties in mind, variational algorithms have been re¬ 
cently developed for particle dynamics [3-5] as well as for 
Vlasov-Maxwell and Vlasov-Poisson systems [6, 7]. 

Action principles in plasma physics (see e.g [8-12] for 
Eulerian formulations of the Vlasov-Maxwell system) are, 
however, still lacking the collisional effects. Emerging in 
the form of a collision operator in the kinetic equation for 
a single-particle distribution function, they describe the 
binary Coulomb interactions between particles and play 
an essential role providing equilibration and increasing 
entropy. Only together, the Vlasov-Maxwell system and 
the collision operator provide a fundamental description 
for plasmas. An action principle for Coulomb collisions 
is thus desired and, in this letter, we provide one. 

The Landau-Fokker-Planck equation, describing the 
collisional evolution of the distribution function in the 
velocity space, is a non-linear integro-differential equa¬ 
tion and, unfortunately, does not exhibit a natural La¬ 
grangian. Thanks to Rosenbluth ]13], the Fokker-Planck 
equation can, however, be formulated as a set of coupled 
partial differential equations: a diffusion-advection equa¬ 
tion for the distribution function / in the velocity space, 
and two Poisson equations for potential functions (/) and 
■0 from which the velocity-space advection and diffusion 
coefficients can be computed. As Ibragimov [14, 15] has 
found a way to obtain Lagrangians and conservation laws 


for any system of partial differential equations of arbi¬ 
trary order, we combine the two ideas and find an action 
principle for Coulomb collisions. 

The letter proceeds by first shortly reviewing the clas¬ 
sical field theory and its extension by Ibragimov for 
our purposes. We then derive a generalization of the 
energy-momentum stress tensor for Lagrangians contain¬ 
ing derivatives up to second order, and apply the theory 
to the Rosenbluth-Fokker-Planck system. We obtain two 
second order Lagrangians and show that the total energy 
and momentum carried by the fields of the extended sys¬ 
tems are conserved. Finally, before summarizing our re¬ 
sults, we discuss the connection of our action principle to 
the Langevin equation for an individual particle. 

As variational integrators have already been con¬ 
structed for specific nonvariational partial differen¬ 
tial equations that admit first order extended La¬ 
grangians ]16], we anticipate the procedure to be applica¬ 
ble to second order Lagrangians, and to the Rosenbluth- 
Fokker-Planck system as well. Finding such a variational 
algorithm is, however, left for a future study. 

Mathematical preliminaries - Throughout this let¬ 
ter, we use the following notation: the independent 
variables (coordinates) in the phase-space-time are de¬ 
noted by a; = ■ ■ ■ i^:™) with m the di¬ 

mension of the phase-space and the phase-space compo¬ 
nents being The dependent variables (fields) 

are denoted by w = {«“}„=! = (m^, . ■., u”) with n the 
number of fields. Partial derivatives of the dependent 
variables with respect to the independent coordinates are 
denoted by the convention it“ = du'^/dx^, which allows 
for a compact notation for a collection of derivatives ac¬ 
cording to U( 1 ) = U( 2 ) = etc. We also 

use the differential 


dxn 


du° 




( 1 ) 


and define the action A[v] for the fields u as the integral 



£(a:,M,U(i),U(2),...)da:, 


( 2 ) 


where C(x, u, M(i), W(2)j ■ • ■) is the Lagrangian (density). 
Summation over repeated indices is assumed. 
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In most of the classical problems in physics that pos¬ 
sess known Lagrangians, the Lagrangian contains deriva¬ 
tives only up to first order, i.e., and the 

Euler-Lagrange equations resulting from the principle of 
stationary action take a familiar form 


dC 

du°‘ 


-D, 


dC 


= 0 . 


( 3 ) 


In this case, the Noether theorem provides conservation 
laws for each symmetry in the Lagrangian according to 


= 0 , 


f = + {r 



( 4 ) 


where the flux is also known as the Noether current. 
The vector fields ^ and rj appearing in have to satisfy 
the invariance condition 


X{C)+CD^{£,^)=Q, (5) 

where the generator for the infinitesimal transformation, 
under which the Lagrangian is invariant, is given by 

^ ik- 

According to Ibragimov [14, 15], one may construct a 
Lagrangian for an arbitrary system of partial differential 
equations of any order (s) as long as the number of de¬ 
pendent fields equals the number of differential equations. 
The recipe is simple: for a set of n differential equations 
for n fields / = {/“}„=! that are defined implicitly by 


and, in case of second order Lagrangians, the flux is 






dC dC 


dC 

du 




( 12 ) 


with IT“ = 77 “ - 

Energy-momentum stress tensor - Before attend¬ 
ing the Rosenbluth-Fokker-Planck system, explicit con¬ 
servation laws for a second order Lagrangian are de¬ 
rived. They are obtained by investigating phase-space- 
time symmetries, and by generalizing the classic defini¬ 
tion of the energy-momentum stress tensor. 

Let us define an infinitesimal transformation with the 
generating vector field X = C^d/dx^ where are con¬ 
stants (constant phase-space-time translations). The flux 
defined in Eq. (12) then becomes where 

the generalized energy-momentum stress tensor^ ap¬ 
plicable to second order Lagrangians has the explicit form 




du\ 


-+uT,Da 


dC 

du' 


fl<T 


dC 
' du 


• ( 13 ) 


ficr 


It is clear that if the Lagrangian is independent of second 
order derivatives, our definition reduces to — 

u'^dC/du'^, i.e., to the standard energy-momentum stress 
tensor for first order Lagrangians. 

We then apply Eq. (11) to find = C‘'C^ 

where the derivative of the Lagrangian with respect to 
coordinate x'^ is understood as an explicit derivative with 
respect to the independent variables. Since the C' are 
arbitrary, the equation holds component-wise, and we 
finally obtain 


Pc{x,f,f(i),...,f(s)) = 0, a = l,...,n (7) 


= = (14) 


the Lagrangian is given by 


£ = 


It is then natural to consider the components of the 
energy-momentum stress tensor separately, and to define 
(8) the energy density 


where /* = {/*“}"=i are the so-called auxiliary fields. 
Further, the adjoint equations, that describe the evolu¬ 
tion of the auxiliary fields, are given by 

i"*(a^,/,r,/(i),/(*i),---) = ^ = o, (9) 

where the Euler-Lagrange operator is defined by 

j Ml■ • -Mj 

In other words, the adjoint equations are found by re¬ 
quiring the action A[f, /*] = / f*°‘Fadx to be stationary 
while computing the variation with respect to the fields /. 

Ibragimov also discusses the generalized conserva¬ 
tion theorem applicable to Lagrangians of arbitrary or¬ 
der. For fields that satisfy the Euler-Lagrange equation, 
6C/Su°‘ = 0, the general theorem can be written as 


( 11 ) 


dC 

e.C-uf— + u?D. 

the energy-density flux 

dC 


dC 

du?„ 


dC 

duZ, 


— u. 


dC 


o dC 
duZ„ ’ 


(15) 


(16) 


the momentum density 

dC 

-PiAD 
^ * du? + * " 


dC 

du?„ 


- u. 


dC 

du?„ 


(17) 


and the stress tensor 
dC 


dc 

du?„ 


— u 


dC 


(18) 


in terms of which, we write equations for the energy and 
momentum explicitly 


Dt£ + D,S^ =Ct, 
DtP^ + =Cj. 


D^f = X{C)+CD^{£,^) 


(19) 

( 20 ) 
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It is clear that a conservation law exists for the energy 
(momentum) related to the fields if the Lagrangian 
does not explicitly depend on time (phase-space coordi¬ 
nates). 

Rosenbluth-Fokker-Planck system - Instead of the 
integro-differential formulation by Landau, Rosenbluth 
found a way to formulate the collisional evolution of the 
distribution function in terms of a coupled system of par¬ 
tial differential equations [13]. Here, for the sake of clar¬ 
ity, we consider only the single species system, which 
is normalized to dimensionless units. Including more 
species is straight-forward. The phase-space coordinates 
X® now denote the three Cartesian velocity coordinates. 

The non-linear Rosenbluth-Fokker-Planck system for 
the distribution function / and for the potentials (j) ^-nd 
ip is then defined by the differential equations 


Ff = dtf - A (/ - AAV® Djf) - (21) 

A = AA</> - /, (22) 

4') = A A^ - (23) 

4') = AAAAV' - (24) 

where we have denoted two different options for the equa¬ 
tion A = 0- We then simply introduce the auxiliary 
fields (/*, AA*) obtain the Lagrangians 

4FP = rA + AA + V®A4. (25) 


As such, contains derivatives of order three and 

derivatives of order four. Addition of total deriva¬ 
tives into the Lagrangian does not, however, alter the 
Euler-Lagrange equations, and we may reduce the order 
of derivatives present in the Lagrangians. 

Neglecting total derivatives, we obtain the second- 
order Lagrangians 

4 fp = Ir {Dt - A<v A) / - i/ (A - AV> A) r 

- - AA AAV® A/ 

- AA A<V - A AAV® - A/ - AV>, (26) 

4 fp = \r (A - A<V A) f-\f (A - AV> A) A 

- AAA<V + A AAA AV - A/ - A/, 

(27) 

and it is straight-forward to verify that both lead to the 
Rosenbluth-Fokker-Planck system when varied with re¬ 
spect to the fields (/*, AA*)- 

Neither £^4 '^m’p depends explicitly on time t 

or velocity x®. It is thus possible to find exact conserva¬ 
tion laws in the field-theory sense simply by computing 


the energy density £1, the energy-density flux S'®, the mo¬ 
mentum density P®, and the stress tensor H®! according 
to their definitions given in equations (15-18). 

Further, since the conservation laws are expressed in 
divergence form and all our fields vanish at the bound¬ 
aries of the infinite velocity space, we immediately obtain 
conservation of the total field-related energy 

A J Sdx^dx^dx^ = 0 , (28) 

and component-wise conservation of the total held- 
related momentum 

A / A dxAxAx'^ = 0, i = I,2,3, (29) 


where the explicit expressions for the energy density and 
momentum density are 

S = £ + i/AA-^AA/, (30) 

A = i/AA - ^AA/. (31) 

Although the dehnitions for energy and momentum 
carried by the extended Rosenbluth-Fokker-Planck sys¬ 
tem are not the physical kinetic energy and momentum 
carried by the distribution function /, the exact conser¬ 
vation laws in our extended continuous system motivate 
finding variational integrators that would provide dis¬ 
crete counterparts for the continuous conservation laws, 
and thus solve the extended system as well as possible. 
As the physical system can be derived from the extended 
one, it is anticipated that the superior properties of vari¬ 
ational integrators, though applied to the extended sys¬ 
tem, would then carry to the physical system as well. 

Adjoint system and Langevin equation - The ad¬ 
joint equations for the fields (/*,</>* A*) derived by 
applying the second order Euler-Lagrange operator 


dL 

du°‘ 


-A, 


dC 

du\ 


aa 


dC 


= F* 


(32) 


with respect to the fields (/, </), "0). For the Lagrangian 
£4 p obtains 


441 ) = - A/* + (AA(A/*) 


+ A (AAAA/*) - A, 

(33) 

4^') = A A A - A - A (/A/*), 

(34) 

441 ) = aAA - AA (a/a/*) > 

(35) 


( 2 ) 

and for the alternative Lagrangian the computation 
gives 

= _aa + (aa(aa) 

+ A(AAAAA)-A-A, (36) 
4’^'^ = A A A - A (/A/*), (37) 

44 ") = A A A A A " ^®^J (A/a A) • (38) 
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Considering then the case of test-particles in a given 
background plasma, the Rosenbluth potentials are pre¬ 
defined and no longer dependent variables: the La- 
grangian becomes simply £r.fp = f*Ff- Noting that 
DjDjip = (j) we define the coefficients 

/r* = -2D,(j) (39) 

(40) 

and write the equations for / and /* as 

Ff = -Dtf - A (/i7) + A A /) , (41) 

f; = Dtf + {D,n + (AA/*) > (42) 

which are the Kolmogorov forward {Ff = 0) and back¬ 
ward (A = 0) equations. The forward equation de¬ 
scribes time evolution of the probability density of a 
stochastic process X that obeys the Langevin equation 
of the Ito type 

dX^ = fi\X)dt + (j^^{X)dW\ (43) 

where dW^ is an infinitesimal change in a Wiener pro¬ 
cess with expectation value E[1T'1] = 0 and variance 
Var)!!^-^] = t. The backward equation, on the other hand, 
describes the evolution of the probability density for t < s 
given end condition f*{x) = p{x). 

Summary - In this letter, we have derived an Eule- 
rian action principle for Coulomb collisions in plasmas 
using the Rosenbluth potentials and the concept of ex¬ 
tended Lagrangians. We also derived the exact conser¬ 
vation laws for the total energy and momentum carried 
by the fields in the extended system. In the case of a 
test-particle population in a given background plasma, 
the action principle was shown to equal the Langevin 
equation for individual particles. 

By providing an action principle for statistical descrip¬ 
tion of Coulomb collisions, this letter, first of all, enables 
for a comprehensive field theory in plasma physics. Sec¬ 
ondly, as the action given is in a form suitable for dis¬ 
cretization, the letter facilitates the development of vari¬ 
ational integrators for the Fokker-Planck collision opera¬ 
tor. Finally, derivation of exact conservation laws allows 
variational algorithms to be benchmarked and validated. 

The author is grateful for the encouraging comments 
from T. Fiilop, J. Candy, I. Pusztai, A. Stahl, O. Embeus, 
G. Papp, A. Bhattacharjee, T. Kurki-Suonio, S. Newton, 
and A. Brizard and for all the fruitful discussions. 
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